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ABSTRACT 

Aims. The dark matter content of early - type galaxies (ETGs) is a hotly debated topic with contrasting results arguing 
in favour of or against the presence of significant dark mass within the effective radius and the change with luminosity 
and mass. To address this question, we investigate here the global mass - to - light ratio T(r) = M(r) / L{r) of a sample 
of 21 lenses observed within the Sloan Lens ACS (SLACS) survey. 

Methods. We follow the usual approach of the galaxy as a two component systems, but we use a phenomenological 
ansatz for T(r), proposed by some of us, that is able to smoothly interpolate between constant M/ L models and a wide 
class of dark matter haloes. The resulting galaxy model is then fitted to the data on the Einstein radius and aperture 
velocity dispersion. 

Results. Our phenomenological model turns out to agree with the data suggesting the presence of massive dark matter 
haloes to explain the lensing and dynamics properties of the SLACS lenses. According to the values of the dark matter 
mass fraction, we argue that the halo may play a significant role in the inner regions probed by the data, but such a 
conclusion strongly depends on the adopted initial mass function of the stellar population. Finally, we find that the 
dark matter mass fraction within R^f / scales with both the total luminosity and stellar mass in such a way that more 
luminous (and hence more massive) galaxies have a greater dark matter content. 
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1. Introduction 

Early- type galaxies (hereafter ETGs) represent the most 
. !^ massive and brightest stellar systems in the universe. 
: Notwithstanding the regularity of their photometric prop- 
5_] ' erties and the existence of remarkable scaling relations, 
' detailed analysis of their structure has been plagued by 
both observational and theoretical shortcomings. Such a 
frustrating situation mainly stems from the lack of a re- 
liable mass tracer able to probe the mass profile outside 
the effective radius. Although the use of planetarey neb- 
ulae have recently improved the observational situation 
allowing the gravitational potential in the outer regions 
dNapolitano et al. 20011 2002, Romanowsky et al. 2003D to 



function with an outer slope s =^ din p/dlnr = —3, there 
is still an open debate concerning the inner scaling of 
PDAiir). The pioneering result pdai ~ of Navarro, 
Frenk and White (1996, 1997) has been questioned by 
other simulations providing either a steeper cusp (e.g. , 
IM oore et al. 1998| [Ghigna et al. 20"00l |Jing fc Suto 2000 



be probed, the well known mass - anisotropy degeneracy 
still represent s a theoretical shortcoming preventing a 
unique reconstruction of the mass profile. 

Numerical simulations of galaxy formation are usually 
invoked as a guide towards understanding the structure 
of the dark matter (DM) haloes. Unfortunately, while 
there is a general consensus on the mass density profile 
PDAiir) being correctly described by a double -power law 
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Fu kushige fc Makino 2001[) or a sha llow profile (e.g 
iPower et al. 20031 INavarro et al. 2004p . Needless to say, 
galaxies are two component systems with the luminous 
matter playing a not negligible role in the inner regions 
which are indeed the most well probed by the data. While 
there is general agreement on the existence and domi- 
nance of DM in the outer gal axy regions, where stars are 
practically absent (see, e.g., Ivan den Bosch et al. 2003a[ 
2003b, 2007), paradoxically the DM content in the inner 
regions is more difficult to interpret notwithstanding 
the availability of a higher number of possible tracers. 
As pointed out by Mamon & Lokas (2005a), the obser- 
vational data claim for a dominant stellar component 
at a radius < 1 Reff- However, the uncertainties on 
which is the stellar initial mass function (IMF), with 
the Salpeter (1995) and Chabrier (2001) as leading but 
not unique candidates, makes it quite difficult to assess 
the general validity of such a result. On large scales, 
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the DM content has been found to be a strong func- 
tion of both luminosity and mass (Benson ct al. 2000', 
IMarinoni .^ Hudson 2002, Ivan den Bosch ct al. 2007) 
with a different behaviour between faint and bright 
systems. Looking for a similar result for the DM content 
within ^ ^Re ff is quite controv e rsial. On one hand , 
some authors ([Gerhard et al. 20Mt IBorriello et al. 2003]) 
argue f or no dependence. In c ontrast, other works 
(e.g., 



2005 Cappellari et al. 2006 



Napolitano et al 
ITortora et al. 2009bp do find that brighter galaxies 
have a larger DM content, while a flattening and a 
possible inversion of this trend for lower mass sys- 
tems, similar to the one observed for late -type galaxies 
( Persic, Salucci fc Ashman 1993P, is still und er analysis 



( Napolitano et al. 20051 ITortora et al. 2009b]) with no 
conclusive result yet obtained. 

Given the controversial situation for local galaxies, it 
is not surprising that the analysis of high redshift galax- 
ies is still more difficult considering that only partial 
information are usually recovered. For instance, rather 
than a full velocity dispersion profile, typically one is 
able to measure only a single value describing the global 
velocity dispersion of the galaxy. This lack of radial 
data extent obviously limits the analysis representing a 
strong obstacle to any attempts to constrain the galaxy 
mass profile. Fortunately, a different mass tracer may be 
available for a subset of ETGs at intermediate redshift. 
Actually, gravitational lensing gives a strong constraint 
on the projected mass within the so called Einstein ra- 
dius (jSchneider et al. 19921 [Fetters et al. 200ip . It is worth 
stressing that such an observational quantity is fully model 
independent although a galaxy model has to be assumed 
to extrapolate the results to other radii. Combining strong 
lensing with the information from galaxy dynamics makes 
it possible to not only increase the number of constraints, 
but also to partially break the mass - anisotropy degener- 
acy ( Treu fc Koopmans 2002a[ [Koopmans fc Treu 2003a[ 
IBolton et al. 2008P thus representing a powerful tool to 
probe mass profiles for lens systems. 

It is worth noting that most of the analyses in litera- 
ture heavily rely on the choice of an analytical expression 
for the density profile of the DM halo. Needless to say, this 
makes the results model dependent with the risk of bias- 
ing them in a difficult to quantify way. Here, in order to 
escape this problem, we follow the approach presented by 
some of us in Tortora et al. (2007) where a phenomena- 
logical ansatz for the global M/L ratio T{r) = M{r)/ L{r) 
has been proposed. Such a model is able to smoothly inter- 
polated between constant M/L models (where light traces 
mass) and the different DM haloes profiles thus providing 
a unifying description of a large class of total mass profiles. 
Fitting the model to the lensing and dynamics observables 
makes it possible to constrain its characterizing parame- 
ters thus giving constraints on the DM mass content and 
its dependence on the luminosity or mass of ETGs at in- 
termediate redshift. As a further step on with respect to 
our original paper, we here improve the analysis by first 
investigating which region of the parameter space is better 
suited to mimic the most used dark halo models. We then 
fit our model to the observed data on lensing and dynam- 
ics in order to individuate, within the plethora of possible 
trends admitted by our parametric choice for M/L ratio, 
the best ones able to match observations. Such an analysis 
then makes it also possible to infer constraints on the DM 



mass fraction and investigate whether any dependence on 
luminosity and stellar mass can be detected in these inter- 
mediate redshift systems. 

The paper is organized as follows. In Sect. [H we intro- 
duce our parametrization of the global M/L ratio and dis- 
cuss the main properties of the effective galaxy model com- 
ing out from the combination of our ansatz for T(r) and 
an approximate analytical deprojection of the Sersic sur- 
face brightness profile. The lensing properties on the model 
are discussed in Sect. [3] where we work out expressions for 
the deflection and Einstein angle. The statistical method- 
ology used to contrast the model against the observations 
and the dataset are presented in Sect. [Sj while a detailed 
discussion of the results is the subject of Sect.[ni We finally 
summarize and conclude in Sect. [7] Along the paper, where 
needed, distances are converted from angular to physical 
units adopting a flat ACDM cosmological model with pa- 
rameters {VlM,^A,h) = (0.3,0.7,0.7) ( jSpergel et al. 2003p . 

2. The model 

The lack of a reliable mass tracer in the outer regions of 
ETGs makes the problem of DM content in these systems 
still open and hotly debated. In the usual approach, a dark 
halo model, described by a given density profile, is added to 
the visible component and then used to fit a given dataset 
(e.g., the velocity dispersion profile or the X-ray emission). 
Needless to say, such an approach is strongly model depen- 
dent and plagued by parameter degeneracies so that the 
constraints on, e.g., the DM mass fraction within the effec- 
tive radius i?e/ / are possibly affected by a theoretical prej- 
udice. In order to overcome such a problem, some of us have 
recently advocated the use of a different strategy. Rather 
than separately modeling the two galaxy components, we 
have proposed to directly constrain its total (luminous plus 
dark) cumulative mass profile M{r). Under the assumption 
of spherical symmetry, it is possible to write it as : 

M{r) = M,{r) + MoMir) = T(r)i(r) 

where Mi, (Mom) is the mass distribution of the luminous 
(dark) component, L(r) the integrated luminosity profile, 
and T(r) is the global M/L ratio. Since L{R) may be re- 
constructed (at least, in principle) by deprojecting the ob- 
served surface brightness, the only unknown ingredient is 
just T(r). On the one hand, should the DM content be neg- 
ligible, which is a borderline case of a constant M/L model, 
one should expect T(r) to remain approximately constant 
to the stellar M/L value. In contrast, T(r) will be an in- 
creasing function of r if the dark halo plays a significant 
role since the total mass M{r) should overcome the stellar 
one in the outer regions. A unified description of these two 
opposite tre nds is provided by th e following phenomenolog- 
ical ansatz ([Tortora et al. 2007^ : 
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with 7/ = r/Reff, 770 = ro/Reff and 
To(l + r;o)'' 



Q+/3 

Vo 



(2) 



While {To,rQ) are scaling quantities, the two parameters 
plays a more significant role setting the asymptotic 
slopes of the global M/L ratio. Indeed, it is : 



T(r) oc 



'^o{r/ro)" r « tq 



(3) 



Eq.([3]) clearly shows that, rather than /3, it is a + 
the important quantity so that, in the following, we 
will parametrize the model by the parameters (a, a + 
(3, rjQ, Te//). It is instructive to look at the following limits : 



lim T(?7) = < 



a>0, a + l3>0 

00 a<0, a + f3>0 

a> 0, a + P <0 

.00 a<0, a + P<0 



(4) 



while for a ~ 0, we get : 



lim T(?7) = To 

r;— ►O 



(5) 



According to Eqs.lg]) and (O, the global M/L diverges if a 
is negative. Considering its definition and that both M(r) 
and L{r) identically vanish in the origin, we shall conclude 
that, for this class of models, M{r) goes to zero faster than 
L(r) which is not expected. Indeed, should visible mass 
dominate in the inner regions, we expect M(r) ~ M^,{r) = 
T* L{r), with the stellar M/L ratio, so that M(r) and 
L{r) decreases with the same rate. In the opposite case, i.e. 
DM is more concentrated than the light, it is possible to 
have M(r) decreasing slower than light so that we cannot 
discard models with negative a. On the other hand, when 
a > 0, we get a null global M/L ratio in the origin which is 
an unphysical result since T(r) can never become smaller 
than the stellar M/L. For this same reason, the case a = Q 
should be accepted provided that Tq > hold. 

One could then argue that the case a > should be 
discarded. It is worth stressing, however, that our ansatz 
in Eq.([T]) is purely phenomenological. We have looked for a 
functional shape that makes it possible to fit as well as pos- 
sible different realistic models suitably adjusting the char- 
acteristic parameters. As for all fitting expressions, we do 
not expect that the fitted approximation recovers the un- 
derlying dark halo model from zero to infinity, but only 
over a limited (although large) radial range. As we will see 
later, this is what indeed happens with our parametriza- 
tion of T(r) notwithstanding the very inner profile shape. 
Motivated by these considerations, we will therefore con- 
sider both a > and a < models. 



As an attempt to constrain (a + (3), we can consider the 
following limits : 



lim T(77) 

ri — ^00 



00 a>0, a + (3>0 

00 a<0, a + /3>0 

a>0, a + /3<0 

lO a <0, a + P <0 



(6) 



which show that models with a + /3 < are clearly unphysi- 
cal. Indeed, in this case, we easily get M(r 00) = 0, that 
is to say the total mass of the system is null. Actually, as we 
will see later, it is not possible to reproduce the T(r) trend 
of any DM model with our ansatz ([1]) when a+/3 < so that 
we will henceforth only consider models with nonnegative 
values of the a + (3 parameter (Tortora et al. 2007[) . In this 
case, the total mass is divergent, but this is a common fea- 
ture for many dark halo models (e.g., the isothermal sphere 
and the NFW profile). In order to avoid this problem, one 
may truncate the halo at the virial radius. 

Further hints on the role of the different parameters at 
work may be obtained considering the logarithmic slope of 
the global M/L ratio. It is: 



dlogT(?7) 
d log 77 

so that : 

dlogT(r/) 
dlog?7 



a + (}ri{ri + rjo) ^ 



(7) 



> 



V>Vc = - [a/ {a + (3)]i]o 



Since we have decided to only consider models with a + 
/3 > 0, we therefore conclude that the global M/L ratio 
is everywhere increasing if a > 0, while the logarithmic 
slope becomes positive only after the threshold radius rjc 
when a is negative. As we will see, both these trends are 
possible for realistic DM models so that we will not set any 
constraint on rj^. It is worth noting that, in a narrow range 
around rj^, the logarithmic slope is approximately null so 
that T(77) is constant as in the mass - traces - light models. 
Should the data at hand probe just this radial range, one 
can reasonably fit them with a constant M/L ratio model. 

2.1. The luminosity density 

The phenomenological ansatz for T(r) is just one side 
of the coin, the other one being the choice of the lu- 
minosity profile. In principle, this quantity may be di- 
rectly reconstructed from the data by first deconvolv- 
ing the observed surface brightness with the PSF of 
the images and then deprojecting it under an assump- 
tion for the intrinsic flattening. Needless to say, such 
an idealized procedure is strongly affected by noise and 
will call for a case - by - case study. Fortunately, ETGs 
are quite regular in their photometric properties. Indeed, 
as well known ((Caon et al. 19931 IGraham fc: CoUess 19971 
|Pru gniel fc Simien 1997), their surface brightness is well 
described by the Sersic profile (jSersic 1968) : 



I{R) = le exp < -bn 



R 



R, 



1/ti 



(8) 
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with R the cylindrical radiusQ on the plane of the sky and 
/e the luminosity intensity at the effective radius Reff- The 
constant 6„ is determined by the condition that the lumi- 
nosity within Reff is half the total luminosity, i.e. 5„ is 
found by solving : 



r(2n,&„) =r(2n)/2 



(9) 



where r(a, z) is the incomplete F function. Although Eq.® 
is straightforward to solve numerically, a very good analyt- 
ical approximation is given by (jCiotti fc Bertin 1999P : 

1 0.009876 

bn ^2n- . 

6 n 

The deprojection of the intensity profile in Eq.® is 
straightforward under the hypothesis of spherical symme- 
try, but, unfortunately, the result turns out to be a some- 
what involved combinations of the unusual Meijer functions 
( [Mazure fc Capelato 2002p . In order to not deal with these 
difficult to handle expression, we prefer to use the model 
proposed by Prugniel and Simien (1997, hereafter PS97) 
whose three dimensional luminosity density reads : 



jir) = jo 



R. 



■eff 



exp 



-bn 



R. 



eff 



l/r, 



(10) 



with 



Jo 



lob^^'-P-^ T{2n) 
2Reff r[n(3-p„)] 



(11) 



Here, Iq = I{R = 0) = /ec''", while the constant p„ is 
chosen so that the projection of Eq. pO]) matches a Sersic 
profile with the sam e values of (n, Reff, l e)- A useful fitting 
formula is given as dMarquez et al. 200 ip : 

0.6097 0.00563 



Pn - 1.0- 

IL n- 

Because of the assumed spherical symmetry, the luminosity 
profile may simply be obtained as : 



L{r) = 47r / r'^j{r')dr' 
Jo 



which, for the PS model, becomes : 



L{r) — Ltx 



7[n(3-p„),6„r;i/") 
T[n{3-pn)] 



where the total luminosity Lt reads : 
Lt = 2T:nb~^"e''"T{2n)IeRlff . 



(12) 



(13) 



^ Note that we have implicitly assumed that the intensity 
I does not depend on the angular coordinates. Actually, the 
isophotes are not concentric circles, but rather ellipses with 
variable ellipticities and position angles so that / = I{R,ip). 
However, in order to be consistent with our assumption of spher- 
ical symmetry of the three dimensional mass profile, we will 
neglect such an effect and, following a common practice, circu- 
larize the intensity profile considering circular isophothes with 
radii equal to the geometric mean of the major and minor axes. 



Note that the total luminosity is the same as the projected 
one for the corresponding Sersic profile as can be immedi- 
ately check computing : 



Lt = 2-k / I{R)RdR 
Jo 



It is worth remembering the definitions of the different spe- 
cial functions entering Eqs.® - : 



r(a,y) 



7(a,y) = / e-'e~'dt , 
Jo 

so that the following useful relations hold : 
r(a) = limr(a,y) = lim j{a,y) , 

y^O y^oo 

l{a,y) = r(a) - r{a,y) . 

As a final remark, let us stress that the total mass of 
the stellar component Mj = Lt may be deter- 
mined from the measurement of the photometric param- 
eters {n, Reff , le) provided that an estimate of the stellar 
M/L ratio is available (for instance, from the relation 
between and the colours or from fitting the galaxy spec- 
trum to stellar population synthesis models). 



2.2. The total mass, density profile and DM fraction 

Combining Eq.ljT]) and ([H]), we trivially get for the total 
mass profile : 



M{r) = ToLt 

= Te//LT 



V 



1 



m/ V Vo/ T[n{3-p„ 

??°(^ + ??o)''7["(3-p«), 
{l + 71of r[n(3-p„)] 



(14) 



Formally, M{r) diverges if a -I- /3 > 0. However, this is not 
a serious problem since, following a common practice, we 
will truncate the model at the virial radius Rmr defined as 
the radius of the sphere within which the mean density is 
'^ciz)j?M{z) with : 

Pm(z) = noMPcritil + zf = {3H^/8TTG)noM{l + zf , 

and, following Bryan & Norman (1998), we approximate 
the virial overdensity as : 



Ac(^) 



187r + 82a;(z) - i^x^iz) 



OM 



where f^oM is the present day matter density parameter 
and, for the fiat ACDM model we are considering, it is : 



x{z) = Qniz) ~ 1 



17om(1 + z)3 
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The total mass is then the virial one obtained by setting 
"n = rivir = Rvir/Reff in Eq.([Tl]) with : 



Rv 



3M„ 



1/3 



and solving with respect to M^r- This gives : 



(15) 



M, 



so that one can conveniently rewrite the mass profile as : 



M(r) = M^,r 

ThT ^ > 

which clearly shows that {Ref fjlei^eff) only works as scal- 
ing parameters through Myi^, while the shape of M{r) is 
controlled by the Sersic index n and the M/L slopes {a, (3). 
Note that the light traces mass case is only achieved by set- 
ting both a = (3 = and "^eff = T*, with T, the stellar 
M/L ratio. Varying these parameters, Eq. (jl6p provides a 
wide range of possible mass profiles thus allowing a rather 
large class of dark halo models to be mimicked. Note also 
that, in principle, one should also check that M{r) is a 
monotonically increasing function of r. To this end, it is 
convenient to consider the logarithmic slope of the mass 
profile reading : 



dliiM 

dim] n7[n(3 - Pn), &n??^/"] 

HV + Vo) + 7["(3 - Pn),bnii^/'' 



+ 



V + Vo 

It is instructive to first consider the case n = A. We get 



(17) 



, dlnM 

1™ ~^^ — 



din A/ 
^™ ~T^ 



a + P 



so that only models with a > and a + (3 > are viable. 
Leaving n as a free parameter, we find : 



dluM 
dhirj 



> 



ario 



Not surprisingly, this is the same constraint we have ob- 
tained when considering the logarithmic slope of the global 
M/L ratio T(r) so that we refer the reader to the above 
discussion. We only stress here that models with a < have 
an unphysical mass profile in the inner regions. As such, the 
behaviour of these models for such low values of r should 
not be considered as a good approximation of the reality. 
Indeed, we get : 



lim M{t], a < 0) 

J7— >0 



±CXD 



which is clearly an unphysical result. However, one can con- 
sider our ansatz for T(r) only as a phenomenological ap- 
proximation working well everywhere but for rj < rjc- Since 



this quantity is typically quite small, hereafter, we will con- 
sider both negative and positive values of a, while only 
models with a + (3 > will be explored to avoid dealing 
with an unphysical vanishing total mass. 

Because of the assumed spherical symmetry, the mass 
density may be straightforwardly derived as : 



p{r) 



1 dM 1 1 dM 



A'Kr'^ dr AttR^j:^ rf drj 
Using Eq. (fTB|) and some lengthy algebra, we finally get : 

M^y2^ 



piv) = 



with : 



A-KRl^^Mvtr 



pi{-n)p2{ii) 



(18) 



Pi = — 



a-3 



Vo 



AA_=(^) fl + ^) 7K3-p„),Mr], (19) 



1 + ^ 
Vo 



exp ( 



(20) 



= (^l + -^)(6„^l/«)n(3-p„) 

7[n(3-p„),6„?7'/"] e^"''''"(21) 



a + 



It is worth stressing that Eq. (|T8|) reduces to the stellar mass 
profile setting a = /3 = and Tej^ = T» as expected, 
while it is positive everywhere for all other values of the 
parameters provided the above constraints on a and a /3 
are set. In order to get further insights on the shape of the 
density profile, one may consider its logarithmic slope : 



s(j]) = 



d In p J] dp 



d In r p dr\ 
Some lengthy algebra finally gives : 

sWi = si(77) -t- 53(77) , 
having defined : 



(22) 



, N g - 3 + (« + /3 - 4) ?7/?7o) , . i/„ 

siW^ T-—, {bn/n)T] ' , (23) 

l + ri/r]o 



S2{v) 



(24) 



X {(3 - pn){l + v/Vt)) + a [1 + (a + /3){ri/ija) /a]} 

+ [n{a + (3)+ am [1 + (a + /3) (77/770) /a] r,'^+^^''} 

X 7["(3-p„),M'/"] , (25) 

V2 = (6„77V")"''"^"ni + '7/^o)+^a7K3-p„),M'/"] 

X [1 + (a + /3)(77 + 77o)/a] exp (6„77i/") . (26) 
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The density profile may be locally approximated as a power 
law with a running slope, i.e. we may write p{ri) ~ rj^^"^^ , It 
is therefore particularly interesting to consider the asymp- 
totic limits of the logarithmic slope. To this aim, we first 
set n = 4 to get : 



lim s(?7, rt = 4) = a — 3 , 

r;— >0 



lim s(ri,n 

rj — >oo 



4) 



-oo 



(27) 



(28) 



so that we get the constraint a < 3 in order that p{t]) is 
not an increasing function of the radius in the inner re- 
gions. Note also that models with a = 3 provide a den- 
sity profile with an inner core thus resembling popular 
cases such as the non singular isothermal sphere. On the 
other hand, in the asymptotic limit i] — > oo, the den- 
sity profile exponentially fades away. Such a behaviour 
should not be surprising being trivially a consequence of 
the exponential decay of the luminosity density profile j{r). 
Although somewhat extreme, such a shape in the outer re- 
gions is not unrea listic. Indeed, recent high resolution N- 
body simulations (jNavarro et al. 20041 IMerritt et al. 2006P 
show that the global mass profile is well approximated 
by an Einasto model (Einasto 19651 lEinasto fc Hand 19891 
ICardone et al. 2005) which indeed predicts such an expo- 
nential decaj0. It is worth stressing that the same results 
are obtained also for other values of n as we have checked 
numerically (for n > 2) so that the constraint on a can be 
used in the following analysis. 

As a final cautionary remark, we stress that the loga- 
rithmic slope 5(77) in Eg. ([22)) refers to the density law of 
the effective model defined by our ansatz for the global 
M/L ratio. As such, 3(77) is not the logarithmic slope of the 
(eventually present) dark halo model. Actually, while in the 
outer regions the stellar density may be neglected so that 
s{r]) may be identified with the soMiv): ™ the inner regions 
things are more involved. In particular, an inner core of the 
effective density profile, i.e. 5(77 — > 0) = 0, does not imply 
that the corresponding dark halo model is a cored one. 

Finally, one of the main outcomes of our analysis is 
the spherical DM content in the core of galaxies, typically 
within 1 Reff- We define the spherical DM fraction within 
the radius r as 



foAiir) = 



MoMir) 



= 1 



M,(r) 



(29) 



MDM{r) + ]VL{r) M{r) 
where Mi,{r) — T^L{r) and MoMif) is DM mass. 

2.3. The surface mass density 

A further quantity to be discussed is the surface mass den- 
sity, which we compute as : 



having defined R = {x^ + y^Y^'^ as usual. It is convenient 
to rewrite the above integral as : 



with r = R/ sinO, so that, inserting Ea. lfTS]) for p{r) and 
defining ^ = R/R^ff, after some algebra, we finally get : 



m - . ""zSZ m 



2Trm]gRl^jjVvtr 



(30) 



with : 



sm( 



X exp 



l/n' 



(31) 



m - 1) 



2TrRl^j n(l + 77o)'3r[n(3-p„)] 



(32) 



Eq. ([30|) clearly shows the role of the different model pa- 
rameters. Indeed, while {Reff,Ie,^eff) only work as scal- 
ing quantities through Eg/ / , the shape of the surface mass 
density profile is governed by the Sersic index n, the slope 
parameters (a,/3) and the transition radius 770- Ideally, one 
can therefore put further constraints on (a, /?) demanding 
that S(^) is a decreasing function in order to be physically 
viable. Unfortunately, such a program is not feasible. As 
Eg. pop shows, there are two terms determining the shape 
of E('^). While the first one scales as it is not possi- 

ble to infer a general scaling rule for E(^) since it has to 
be evaluated numerically. We have therefore investigated 
the shape of E(if) by evaluating it numerically for differ- 
ent values of {n, a,a + P, 770). It turns out that the surface 
density is a monotonically decreasing function for ^ > 0.01 
whatever are the values of the model parameters. In the 
very inner region, however, it is sometimes possible to get 
a cored profile or an unphysical increasing E(^). Since we 
are unable to work out any analytical rule to identify these 
models, we prefer to not put any further constraints on the 
model parameters. Note, however, that, also in this dis- 
turbing circumstances, S(^) is bad behaved only in a very 
inner region where our phenomenological ansatz may break 
down. Nevertheless, as we will see later, this has no effect 
on the estimate of the quantities of interest. 



/"OO 

E(i?) = 2 / p{r)dz 
Jq 

^ Actually, these are DM only simulations so that the addi- 
tion of a baryonic component changes the global mass profile. 
Nevertheless, in the very outer regions, the stellar component 
may be neglected so that it is likely that the global mass profile 
behaves as the DM one thus recovering the exponential decay 
predicted by these simulations. 



2.4. Projected mass and DM fraction 

A quite interesting derived quantity, which will be useful 
when contrasting models against observations, is the pro- 
jected mass within a circular aperture of radius R given 
by: 

MprojiR) = 2tt [ i:{R')R'dR' = 2TTRljf ^i^')Cd^' 
Jo Jo 
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which, for our model, reduces to : 



pro J 



(0 



nT[n{3 - Pn)] \ 1 + ^70 



.(33) '^'^p 





AttG 



3-^2 (^ap) 



rj{r)M{r)dr 



For the Prugniel- Simien model we are adopting for the 
luminosity component, the projected mass is given by: 



^ r(2n,b„^i/") 
r(2n) 



so that the projected DM mass fraction reads : 



(34) 



fproj {£,) — 1 



= 1 - 



1 



nr[7i(3 -p„)] 
l-r(2n,b„C'/")/r(2n) 



(35) 



It is worth noting that the projected mass fraction is typi- 
cally higher than the spherical mass fraction foMif)- This 
is a consequence of the integration along the line of sight 
which is performed in order to evaluate Mproj{R)- Indeed, 
the stellar component contributes mainly to the inner re- 
gion, while the DM halo extends over a larger part of the 
integration domain. As a result, fprojir) is higher than 
foAiir), with the ratio fDM{r)/fproj{r) depending both 
on r and the concentration of both the stars and dark halo 
density profiles. 



3. Dynamical and gravitational lensing properties 

While photometric observations probe the light distribu- 
tion (hence giving constraints on the luminosity density), 
dynamical observables (e.g., velocity dispersion) and grav- 
itational lensing are excellent probes of the mass profile. 
In the framework of our phenomenological scenario, we can 
tell that dynamical and lensing observables make it possi- 
ble to constrain the global M/L ratio T(r). Therefore, it is 
interesting to derive these main properties for our model. 



3. 1 . Velocity dispersion 

A widely used probe to constrain the model parameters is 
represented by the line of sight velocity dispersion luminos- 
ity weighted within a circular aperture of radius Rap- This 
can be easily evaluated as : 

^2 _ j^-'' I{R)al,{R)RdR 
j""" I{R)RdR 

with uios{R) the velocity dispersion projected along the line 
of sight. This latter may be found solving the Jeans equa- 
tion provided an expression for the anisotropy parameter 
is chosen. In order to not further increase the number of 
unknown quantities in the problem, we assume isotropy in 
the velocity space. In this case, it is possible to arrive at 
the following expression for dap (jMamon fc Lokas 2005aP : 



with 

L2{R) 



(r2 - Rl\^l^j{r)M{T)r~^dT 



I{R')R'dR' 



(36) 



(37) 



the projected luminosity profile. For the model we are con- 
sidering, after some algebra, we finally get : 



„2 

"ap.th 



GTqLt 6"^^"''"^r(2n) 



^ap{£,ap: P) 



GT 



eff^T 



Reff 3nT[n{3 - pri)] 

1ap(£,api P) 



'7o 



(38) 



where we set ^ap = Rap/ Reff , being Rap = 1.5 arcsec the 
aperture radius (equal for all lenses because (Tap is measured 
with the same instrumental setup) and have defined : 



7-(i) 

-^ap 



t(2) 
-^ap 



(39) 



t(i) 

-^ap 



IQ 

X exp I 



'-P-drj , 



1 - 



r[n(3~p„),M'/" 
r[n(3 -pn)] 



t(2) - 

-^ap 



1 - 



^ 1 + ^ 

Cap \'7o/ V ^70 
X exp(-M'/")(^'-0'/'^ 



(40) 



r[n(3-p„),M'/'^ 
r[n(3 -pn)] 

(2+P"'d?7 . (41) 



For given values of the model parameters, Eq. lfSS]) allows 
us to predict the value of (Top. However, because of the 
uncertaintits on both the photometric parameters and the 
stellar M/L ratio, the predicted CTap will be affected by an 
uncertainty which we compute by the simple propagation of 
errors. We do not report here the final expression for sake 
of shortness, but caution the reader that the uncertainty 
on i?e/ / must be propagated also in the denominator term 

7(2n, 6„^op"), while we neglect the contribution to the in- 
tegral for simplicity. 



3.2. Lensing observables 

Gravitational lensing is a powerful tool to further constrain 
the space of parameters. After having determined the sur- 
face mass density in the previous section, here, we define the 
deflection angle, that allows a model independent estimate 
of projected mass at Einstein radius to be determined. 



8 



V.F. Cardone et al.: Global M/L of SLAGS lenses 



3.2.1. The deflection angle 

A key role in the determination of the lensing properties of 
a given model is played by the deflection angle. For a spher- 
ically symmetric lens, this reads ([Schneider et al. 1992P : 



R 



R'dR' 



where the critical surface density Serif is : 



with Ds, Dl, Dls the observer- source, observer - lens and 
lens - source comoving angular diameter distances, respec- 
tively. Inserting Eq. ([30|l , the deflection angle for our model 
is easily computed as : 



= 1) Se«t ^ 

e 5(e = 1) 

where we have definecO : 
Jo 



(42) 



(43) 



2R, 



eff 



^eff 5(C=1) 



2R 



eff 



TqLt a(C = 1) 



2R. 



eff 



27rT,c„tRlff m]^T[n{3 - pn)] 



2^Tl:cr^tRlff n(l + r;o)'3r[n(3 - p„)] 



(44) 



Looking at Egs. iP^ and it is immediately clear that 
the Sersic index n and the model parameters (a,/?, 770) de- 
termine the shape of the deflection angle, while the effective 
radius i?e/ / , the total luminosity Lt and the mass param- 
eter Te// only work as scaling quantities. 

3.2.2. The Einstein angle 

Should the source be a distant pointlike object (such as 
a quasar), one may observe multiple images and then use 
their positions to constrain the model parameters. On the 
other hand, when the source is an extended object approx- 
imately aligned with the lens along the line of sight, the 
formation of arcs takes place. The angular radius of the 
circle individuated by these arcs is the Einstein angle Re- 
Referring the interested reader to the literature for more de- 
tails, we here only remember that the Einstein angle may 
be determined by solving : 



aiis) = Re = iEReff, 



(45) 



A dimensional remark is in order here. The deflection angle 
is expressed in arcsec so that, although the effective radius -Re// 
is a length, in the right hand side of Eq. (I42|l , it must be expressed 
in angular units. 



witllfl = RE/Reffj so that a measurement oi Re (in 
arcsec) provides a constraint on the model parameters 
through the deflection angle d(^£;). Moreover, comparing 
the definition of MprojiS,) and q;(^), it is immediate to get : 



MprojiRE) 



— irR^Y^crit 



<eR, 



eff 



(46) 



with Reff now expressed in linear units. Ea. (|46p allows 
then to estimate the projected mass within the Einstein 
radius provided that T,crit is known, i.e., the lens and source 
redshift have been measured. 



4. Testing the model global M/L ratio 

Our proposed parametrization for T(r) was motivated by 
the need to have a versatile expression able to mimic a 
large set of galaxy models. Since we assume that the stellar 
component is fixed and described by a PS density profile, 
mimicking a large set of models is the same as mimicking 
the global AI/L ratio predicted by different dark halo mod- 
els. Needless to say, such an ambitious goal is difficult to 
reach so that we expect that Eq.((T]) does a good job over 
a limited radial range. It is therefore worth exploring how 
large is this range and how good and versatile is our ap- 
proximated parametrization. 

To this aim, it is convenient to resort to the following 
general expression for the mass profile of the dark halo : 

MDM{r) = M^ir X ^J.{r/rs) = M^ir x n{ri/ris) 

with Mmr the virial mass and Ts — Reffijs a- scalelength. 
As a function of the adimensional length 77 the DM mass 
fraction then reads : 



MDAliv) 



Normalizing with respect to fe 
using Eq. p2|l , we finally get : 



1 - fe fJ-i'^/Vs) 



foKiiR = Reff) and 



/dm (77) = 1 + 



fe KV) 



1 + r[n(3 - p„), 6„7yi/"]/r[n(3 - p„)] 



1 + T[n{3 - pn),bn]/T[n{3 - pn)] 
which may then be related to the global M/L ratio as : 



(47) 



T(r;) 



1 - foAliv) 



(48) 



having assumed a constant stellar M/L ratio. Ea. (|T7)l 
makes it possible to compute the DM mass fraction and 
hence the global M/L ratio through Eq. (|48|) provided that 
an expression for fi{r]/r]s) is given. 

To test the validity of our phenomenological ansatz, 
we start considering the most popular galaxy profiles, 
i .e. NEW (Navarro et al. 1996, 1997), the Einasto 
Pinasto 19651 lEinasto fc Haud 19891 INavarro et al. 2004[ 
ICardone et al. 2005|) . the nonsingular isothermal sphere 



* Here, we recast the definition for deflection angle and pro- 
jected mass introducing the adimensional radius ^e. 
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and the Burkert (1995) models. Using the formalism in- 
troduced above, the cumulative mass profile for each case 
is assigned by the following expressions : 

fJ-NEwiv/Vs) = log (1 + 77/77^) - (?7/?7s)(l + ?7/77s)"\ 

r[2ncM,rf„(r//r?s)'/"""] 
T{2nDM) 

fJ-Nisiv) ^V/Vs - arctan (77/775) , 

In [1 + iv/Vs)^] , In [1 + iv/Vs)] - arctan (ij/tJs) 
MB«r(?7) = + ^ 

We parameterize NFW profile assigning first the virial 
mass Mj;ir and then estimating the scalelength as — 
Ryir/c using the following numerically motivated relation 
(jBullock et al. 20011 [Napolitano et al. 2005[ ) to set the con- 
centration parameter : 



16.7 

l + z 



10"/i-i M, 



-0.125 



also taking care of the large scatter. It is worth noting 
that, by virtue of the c-M^r relation, the NFW profile 
reduces to a one - parameter model. In contrast, since no 
numerically or observationally motivated relation between 
("n-DM, 'l^s^ Myir) is known, the Einasto model is assigned by 
three parameters that we choose in such a way to have a 
realistic density profile. A similar discussion also holds for 
both the NIS and Burkert models which are, however, two- 
parameters profiles. As a further remark, we stress that, 
while the NIS and Burkert models possess an inner core, 
the NFW and Einasto models are cuspy even if with differ- 
ent logarithmic slope for 77 << 77^. Moreover, the outer slope 
of the density profile runs from s = —2 for the NIS model, 
to s = —3 for both NFW and Burkert cases, to s = —00 
(i.e. exponential decay) for the Einasto parametrization. 
As such, by considering these four models, we are spanning 
a wide range of possibilities thus allowing us to test how 
good our phenomenological proposal for the global M/L 
ratio T(r) mimics a large class of different profiles. 

Once the halo parameters have been set, we add a stellar 
component fixing the photometric parameters (ti, i?e//, ^e), 
the stellar M/L ratio T* and the redshift z so that they 
are equal to those of a randomly chosen lens. We then fit 
Eq.([T]) to the resulting T(7;) profile and repeat this exer- 
cise for 1000 random realizations of each halo model. The 
instructive results of this investigation are schematically 
summarized below. 

- Eq.([T]) fits the T(r) of the tested models with a 
very good precision and over a large radial range, 
namely 0.01 < 77 < 10. Denoting with T fit our 
proposed function and defining A = |1 — T/ii/T|, 
we get that the maximal deviation is lower than 15% 
(depending on the model). In particular, for NFW is 
0.06% < Ar„,s < 1.8% and 2% < A„iax < 8%. The 
Einasto model is fitted with Arms ^ 0.5%, while the 
worst performances are obtained for the Burkert one for 
which we, however, find a still satisfactory Arms 6%. 



The fitted model parameters (a, /3, 770, Te///T*) depend 
on the halo parameters so that a general rule cannot be 
given. However, we note that quite small and negative 
values of a are clearly preferred for the NFW-I-PS 
input model. This is somewhat surprising since this 
would give rise to models with an inner asymptotic 
slope of the global density sq ~ a — 5 < —3, contrary 
to popular models having sq — —1. However, as stated 
above, in the very inner regions our model does not 
fit anymore the NFW -|- PS one so that the results on 
So can not be trusted upon. For this same reason, we 
do not care about obtaining still more negative values 
of a. Indeed, for the other models, a < cases are 
still clearly preferred. As a conservative estimate, we 
therefore conclude that, in order our parametrization 
fits well different dark halo models, one must have 
negative a, with values as low as a ~ —0.40. 

The asymptotic outer slope a + P typically takes values 
of order 1-2, even if one cannot exclude values as 
high as a + P ~ 5. For example, if considering the 
NFW + PS, a -h /3 is strongly peaked on the ^ 2, 
and the distributions is markedly asymmetric with 
long tails towards larger values. In contrast, a + (3 
has a similar distribution, but shifted to the smaller 
a + P 0.5 peak, for cored models. As a general rule, 
however, a + P < never occurs in agreement with 
our above discussion showing that such models have an 
unphysically decreasing mass profile. 

The scaling quantities Te///T* and log 770 are in the 
range 1.0 < Te///T* < 2.5 and —0.5 < logT^o < 1.5. 
In particular, if considering the NFW -I- PS, Te///T* 
and log 770 are strongly peaked on ^ 1.5 and 0.5 
respectively. Similar values are obtained for the other 
dark halo models considered, a result which is not fully 
unexpected. Indeed, Te///T* and log 770 are mostly 
related respectively to the virial mass Myir and DM 
mass fration /e and to the scaled virial radii Ryir/Reff- 
Since these quantities take similar values for all models, 
we indeed expect to find a not too large scatter among 
different models which is indeed what we find. 

Changing the c- M^r relation or its scatter does not 
alter the result on the validity of our approximation for 
NFW -I- PS profile, i.e. Arms and Amax span the same 
range as above for 77 covering the same radial range. 
However, the distribution of the fitted parameters 
may change significantly. In particular, higher values 
of a may be obtained although negative a are still 
preferreqj. Considering that similar values for a can be 
obtained by fitting T fn to models other than NFW, 
we argue that it is not possible to resort to a to 
discriminate between, e.g., the NFW and the Einasto 
models, or cusped and cored profiles, unless one has a 
precise determination of the c- M^r relation. 



^ Note that the relation between the concentration c and a is 
difficult to extract. Indeed, since we set a as a result of fitting 
our phenomenological model to an input stellar plus dark halo 
profile, its value also depends on the fitted values of the other 
parameters, (a + l3,r]o). As a result, we are therefore unable to 
say whether a larger a implies a higher or smaller concentration. 
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These results clearly show that our proposed parametriza- 
tion of the global M/L ratio indeed makes very well its job. 
It is indeed able to reproduce with very good accuracy the 
expected behaviour of M{r)/L{r) whatever is the dark halo 
density profile over a radial range spanning three orders of 
magnitude. As a final cautionary test, we check that using 
this very good approximation does not introduce any bias 
in the lensing deflection angle. As an example, we can con- 
sider the case of the NFW model and comparing the value 
of a{£) computed using the analytical NFW deflection an- 
gle and the numerically evaluated Eq. ([l^ . We repeat this 
test for 1000 random realizations and estimate the ratio 
OiNFwiO/otfAO for ^ = 1 and £, = S,e, using for the 
observational value for the lens used to set the stellar pa- 
rameters. It turns out that both these ratios deviate from 
unity less than 1% confirming our expectation that no bias 
is induced by the use of our ansatz. 



5. Model vs observations 

The phenomenological ansatz for T(r) we are proposing is 
characterized by four parameters, namely the two asymp- 
totic sloped, the logarithm of the scaling length log 779 and 
the global M/L ratio at the effective radius ^eff- To these 
four quantities, we have to add the stellar M/L ratio T* 
thus leading to five the number of unknown quantities to be 
constrained. Needless to say, tackling such an issue is a quite 
daunting task. A possible way out could be contrasting the 
model with kinematic observations, such as the velocity dis- 
persion profile. To this aim, the data should cover a radial 
range wide enough to sample with sufficient detail both the 
inner and the outer regions in order to constrain the two 
asymptotic slopes. Moreover, the measurement errors and 
the sampling should be extremely good to overcome the 
problem of parameters degeneracies in a five dimensional 
space. All these observational requirements are quite de- 
manding so as to be satisfied only by a handful number 
of nearby galaxies. Actually, we are interested here in lens 
galaxies with typical redshifts of order 0.1 — 0.5 so that 
measuring their velocity dispersion profile with the above 
precision and sampling is an unrealistic task. Nevertheless, 
for each galaxy, we still have two observable quantities that 
can be used, namely the Einstein angle Re and Gap, the 
line of sight velocity dispersion luminosity weighted in a 
circular aperture of radius Rap- Needless to say, it is im- 
possible to constraint a five dimensional parameter space 
with only two datapoints so that a case - by - case analysis 
is not possible. In order to overcome such a problem, we 
have therefore to reduce the number of unknowns and in- 
crease the number of observed datapoints. To this aim, we 
therefore first get an estimate of the stellar M/L ratio T*. 



5.1. Stellar M/L 

We start assembling a library of synthetic stellar pop- 
ulation models obtained through the Galaxev code 



® Hereafter, we will consider a + (3 rather than /3 as model 
parameter since it is this former quantity that sets the asymp- 
totic slope of the global M/L ratio. Moreover, we have physically 
motivated constraints on a and a+/3 so that the ranges to be ex- 
plored for these parameters is immediately set. We also change 
from 7^0 to log?7o as scalength parameter in order to explore a 
wider range for this quantity. 



(|Bruzual &: Chariot 20031) varying the age of the popula- 
tion, its metallicity and time lag of the exponential star 
formation rate and assuming a Chabrier (2001) initial mass 
function (IMF). Then, we use the tabulated (w, g, r, i, z) ap- 
parent magnitudes (corrected for extinction) of each lens 
to fit the above library of spectra (suitably redshifted to 
lens redshift) to the colours, thus getting the estimates re- 
ported in Table 1. Note that these values may be easily 
scaled to a Salpeter (1955) or Kroupa (2001) IMF by mul- 
tiplying by 1.8 or 1.125 respectively so that we can explore 
other IMF choicefl In order to get the uncertainties, 
we use a Monte Carlo - like procedure generating a set of 
colours from a Gaussian distribution centred on each mean 
colour and standard deviation equal to the colour uncer- 
tainty. Fitting the colours thus obtained to the synthetic 
spectra for each realization, we generate a distribution of 
fitted parameters. The median and median scatter of such 
a distribution is finally taken as an estimate of T* and its 
uncertainty (see Tortora et al. 2009 for further details). 

5.2. Resorting to universal parameters 

Having reduced by one the number of parameters, we now 
look for a way to increase the number of constraints. To 
this aim, we could stack all the lenses in a single sample. 
However, while one can reasonably argue that the same 
functional expression for T(r) describes the global M/L ra- 
tio of all the lenses, its characterizing parameters are likely 
to change on a case -by -case basis. Indeed, even if we as- 
sume that the DM halo has a universal profile, the details 
of the baryonic assembly may lead to different model pa- 
rameters. In order to partially alleviate this problem, we 
therefore reparametrize T(r) in Eq.([T]) in terms of quanti- 
ties that are likely to be (at least, to first order) universal. 

To this aim, it is worth remembering that previous 
works on fitting galaxy models to the lenses using the con- 
straints from the Einstein ring and the velocity dispersion 
concordantly suggest that the total mass profile at R = Re 
is well approximated by a singular isothermal sphere (see, 
e.g., |Treu fc Koopmans 2004 and refs. therein). Actually, 
rather than constraining the global mass profile, such an 
analysis essentially probe the shape of the density pro- 
file only at the Einstein radius so that one can argue that 
the logarithmic slope of the density profile at Rs/Reff is 
the same for all lenses. Motivated by these previous liter- 
ature results, we therefore assume that se = s(Re / Ref f) 
is the same for all lenses and use this quantity as model 
parameter instead of log 770- To this aim, for given values of 
{a,a + 13, Se), we numerically solve : 

S{RE/Reff) = SE 

with respect to log r/o using Eqs. iP^ - (^5)) to compute 3(77). 
We then note that the DM mass fraction at the virial radius 
is likely to be a universal quantity. Indeed, considering that : 



1 



MJJI^ 



1 



T* Lq 

M^ijjj. 



and setting M^{R^ir)/Myir = esFfb with /b = ^It/^M 
the cosmic baryon fraction ( Spergel et al. 200"3j ) and esf 



^ While this is correct for a Salpeter IMF, since it differs from 
a Ghabrier IMF only for the low mass slope and predict very 
similar colours, for a Kroupa IMF this is not strictly true, but 
we could assume the scale factor above as a good approximation. 
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the star formation efficiency, we can assume that, notwith- 
standing the details of the star formation process, esp is 
roughly the same for all the lenses in the sample. Although 
not strictly true, this is, however, a resonable approxi- 
mation since all the lenses we are considering are ellip- 
tical galaxies spanning a limited range in both redshift 
and luminosity. As such, we will hereafter assume that 
= Myir/T^LT is a Universal quantity to set the T^ff 
parameter as : 



-eff 



r[n(3 -Pn)] 



(49) 



with Myir = fiyirT^^Lx- Note that, here, we are assum- 
ing that M^{Ryir) is the total stellar mass T^,Lt which, 
although not strictly true, is a very good approximation. 

Unfortunately, we are unable to find other lens related 
parameters that can be considered universal so that we will 
hereafter parameterize our T(r) ansatz with the slope pa- 
rameters a and a + P, the logarithmic slope at the scaled 
Eintcin radius se — s{Re/ Ref f), and the mass ratio at 
the virial radius = /T *Lt assuming that these 
are universal quantities. In order to explore the impact of 
this theoretical hypothesis, we will fit our model to the full 
lens sample and to four subsamples containing respectively 
five, six, six, four lenses binned according to the absolute 
V magnitude My - Comparing the constraints on the model 
parameters obtained by the different fits makes it possible 
to look for an eventual dependence on the lens luminosity 
thus giving an a posteriori check of our a priori assumption. 

5.3. Statistical analysis 

As shown by Ea. (|15)l . a measurement of ^e may be con- 
sidered as a measurement of the deflection angle a{S,E) so 
that stacking together many galaxies with different values 
of ^E makes it possible to reconstruct the deflection angle 
profile. Similarly, since Uop is weighted within an aperture 
of fixed angular radius Rap but different normalized radii 
£,ap = Rap/Reff, this is the same as tracing the luminos- 
ity weighted line of sight dispersion profile. The agreement 
within the model and the data may then be optimized by 
maximizing the likelihood function : 



£(p) oc exp 



xLs(P) +Xdyn(P) 



(50) 



with p = {a,a + P, se, ^J'v^r) the set of model parameters 
and Xiens iXdyn) ^ merit function for the lensing (dynam- 
ics) related observables. Maximizing C{p) gives us the best 
fit model parameters, while the isolikelihood contours pro- 
vide constraints on the parameter space that are difficult 
to quantify because of the four dimensional nature of this 
space. In order to carry out the maximization of £(p), we 
implement a Markov Chain Monte Carlo algorithm which 
efficiently samples the regions of high likelihood. The points 
in the chain may be also used to extract constraints on 
a given parameter pi marginalizing over the other ones. 
To this aim, we consider the histogram of the pi values 
along the chain and, following a Bayesian approach, take 
the median of this histogram as the best estimate of the pa- 
rameter Pi, while 68 and 95% quantiles of this distribution 



give us the 68 and 95% confidence ranges for the parame- 
ter Pi. Note that, since the distributions may also be not 
Gaussian, the median and the mean value may be different, 
but we still retain the median as best estimator in agree- 
ment with the Bayesian framework. As a further remark, 
let us stress that the best fit parameters pbf may also sig- 
nificantly differ from the maximum likelihood values Pml 
because of parameter degeneracies. Actually, in order to 
have Pbf — Pml, ^{p) should be written as the product 
of the marginalized likelihood functions. Moreover, these 
latter should be Gaussian so that they are centred on the 
maximum likelihood values (therefore coincident with the 
median values). Actually, both these conditions are rarely 
met so that having pbf Pml is quite common. 

The Markov Chains may also be used to extract con- 
straints on some interesting derived quantities. To this aim, 
let us consider a generic function y — g{p). Evaluating y 
along the chains, we can build up the histogram of its values 
and then use the median and quantiles of this latter to infer 
the Bayesian confidence levels on y. Such a procedure will 
be used for the scaling model parameters logr/o and T^ff 
and for the projected and spherical DM mass fractions for 
each lens in the sample. 

5.3.1. The lensing merit function 

Dealing with the Einstein angle as a constraint is not 
straightforward. Indeed, for a given set of model param- 
eters, £,E should be find solving Eq. (|15)) so that a lot of 
computing time should be spent to explore the four di- 
mensional parameter space. Moreover, the photometric pa- 
rameters (n, Ie,Reff) entering Eq. ([i5)) are affected by their 
own uncertainties so that one should solve this relation 
many times by, for instance, bootstrapping the values of 
these quantities and looking for the corresponding uncer- 
tainty induced on the predicted £,e- Fortunately, we can skip 
all these complications by resorting to the projected mass 
within the Einstein radius whose theoretical and observed 
values may be estimated from Eqs. ([55)) and pS)) . 

Because of the measurement errors on 1^ and Ref / and 
the uncertainty on Tf,ff coming from the one on T*, the 
predicted Mp'^^j{£^E) is affected by an uncertainty which we 
naively quantity using error propagation as : 



<5M^™,(^b) 


A/*'* ■ 

proj 


fe) 


with : 




5Lt 


In 10 











-eff 



(51) 



lnl0./(251ogi?e//)2 + ((51og/e)2 



SAL, 




Note that we have here neglected the contribution of the er- 
ror on Re entering through the integral in Eq. ([55)) . To take 
this into account, we should compute the numerical deriva- 
tive (iM*J?Q (^)/(if for each choice of the model parameters. 
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Table 1. Photometric and lensing observables and estimated stellar M/L ratio for the 21 SLACS lenses. 



Such a computationally expensive procedure may be how- 
ever skipped since, being Re measured with a quite low 
error, the integral function may be well approximated as a 
constant over such a narrow range. As a consequence, the 
corresponding error term is safely negligible with respect to 
the ones due to the measurement uncertainties on /g and 
Reff- We stress, however, that neglecting this term makes 
our analysis more conservative since we are actually assum- 
ing that the final uncertainty on M'^^^IS^e) is lower than its 
true value. The error on Re is instead taken into account 
to compute the one on the observed projected mass being : 

<5M;,%.(^B)/M;,%(Ci;) = 25Re/Re . (52) 
We finally define the lensing merit function as : 



AT 

xLs(p) = 

1=1 



proj. 



eie 



(53) 



where the total uncertainty e;ens,i is given by the sum in 
quadrature of the theoretical and observational uncertain- 
ties estimated from Eas. ((5T|) and ((52|) and the sum is over 
the M lenses in the sample. 

5.3.2. The dynamics figure of merit 

Denoting with cjap.pbs the observed value, we define the fol- 
lowing figure of merit for the dynamical data : 



Xdyn(P) 



E 

i=l 



a 



ap,obs 



,(0 

ap,th 



(p) 



£dyn,i 



(54) 



where e^^^n is obtained by summing in quadrature the the- 
oretical and measurement uncertainties, while the sum is, 
as usual, over the TV objects in the sample. Note that, al- 
though Rap is the same for all galaxies, the different mea- 
surements of (Tap actually probe a range of scaled radii ^ap 



since the effective radius changes on a case -by -case basis. 
As a final remark, it is worth stressing that cTap.th oc '^^f^f, 
while Mp'^^j{^E) oc Tg//. It is this diverse scaling with T^ff 
(and the other model parameters too) which makes com- 
bining lensing and dynamical data an efficient tool to par- 
tially break the degeneracies among model parameters thus 
strengthening the constraints we can obtain. 



5.4. The data 

The sums in Eas. l|55|) and run over the TV objects in 
the sample. It is therefore vital to specify what this sample 
is made of. We consider here the 21 lens ETGs reported in 
Gavazzi et al. (2007) whose main properties are summarized 
in Table 1 where columns are as follows: 1. name (without 
the prefix SDSS); 2. lens redshift; 3. source redshift; 4. ab- 
solute V magnitude (incremented by 51og/i); 5. effective ra- 
dius (in arcsec); 6. Einstein radius (in arcsec); 7. aperture 
velocity dispersion (in km/s); 8. ratio between the aperture 
and the effective radii; 9. estimated V-band stellar M/L 
ratio with its error. Note that we do not have colours for 
lens J223840. 2+075456 so that T* cannot be evaluated and 
this object will be excluded by the analysis. These lenses 
have been observed in the framework of the SLACS sur- 
vey which aims at confirming through ACS imaging can- 
didate lenses spectroscopically identified within the SDSS 
catalog (see Bohon et al. 2004, 2006, 2008 for further de- 
tails). The SDSS velocity dispersions have been measured 
within a circular aperture of fixed radius Rap = 1.5 arcsec 
and spans the range 196 < aap < 333 km/s with a mean 
square velocity (cr^p)^/^ ~ 248 km/s. The lens galaxies have 
a mean redshift (zl) ~ 0.22, but it is worth noting that 
the redshift range covered (0.063 < zl < 0.513) probes 
almost one order of magnitude (even if quite sparsely). 
The measured surface brightness has been fitted by a de 
Vaucouleurs (1948) model so that n will be fixed to 4 in the 
following. Typical uncertainties on Reff and /g are quite 
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small so that, following Bolton et al. (2008), we will set 
(^log/e,(51ogi?e//) = (0.020,0.015) for all galaxies in the 
sample. Gavazzi et al. (2007) provides also absolute V band 
magnitude corrected for filter transformation and Galactic 
extinction and K and evolution corrected to a common red- 
shift z = 0.2. When fitting Einstein radius and velocity dis- 
persions, however, we need the luminosity at the lens red- 
shift. This can be easily estimated as (Gav azzi et al. 2007[) : 

logLy(zL) = \ogLv{z = 0.2) - 0.40(zL - 0.2) 

with : 



logLy(z = 0.2) = dex - 



Mv - Ml 



V,0 



2.5 



having set dex(a;) — 10^ and denoted by My and Mv,q — 
4.83 the absolute magnitudes of the lens at z = 0.2 and 
of the Sun, respectively. We further estimate the effective 
intensity as (Gra ham fc Driver 2005P : 



^6^^"r(2n) 



with the effective surface brightness given by : 
{h)^Lv(zL)l2^Rlff . 



(55) 



(56) 



While Gap is directly measured from the lens spectrum, the 
Einstein radii are determined by a parameterized proce- 
dure. The lens potential is modeled as a singular isothermal 
ellipse, while a non parametric technique is employed to re- 
construct the source surface brightness profile. The Einstein 
radii is then estimated so that the lensed source matches 
the observed arc -like features. Typical uncertainties in 
are of order 0.05 arcsec with a small variation among the 
lenses. Considering the different Re values, we naively set 
6Re/Re = 0.05 for all the sample objects. Table 1 summa- 
rizes the relevant quantities for our analysis. 

6. Results 

The sample of 21 SLAGS lenses represents the input dataset 
needed to constrain the four model parameters (a, a + 
(3 , s E , ^J'vir) through the Bayesian likelihood analysis de- 
scribed above. As a first test, motivated by the above dis- 
cussion, we fit the model to the full sample without any 
binning in luminosity thus giving us a quite large set of 
constraints (namely, 2x21 = 42 observable quantities vs 4 
parameters). Nevertheless, some care is needed when exam- 
ining the results of the Markov Chain analysis. 

In order the results to be reliable, one should carefully 
check that the chains have reached convergence, that is to 
say that the chains have fully explored the regions of high 
likelihood. Should this be the case, one should see the points 
of the chains for each single parameter oscillate around an 
average value or, put another way, the histograms of the 
parameters be single peaked (eventually, with a Gaussian 
shape). In order to test the convergence of the chains we 
resort to the test described in Dunkley et al. (2005, see also 
Dunkley et al. 2009) which is based on the analysis of the 
chain power spectra. Before checking for convegence, how- 
ever, we first cut out the initial 30% of the chain in order to 
avoid the burn in period. Moreover, to reduce spurious cor- 
relations among parameters induced by the nature of the 
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mean 
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95% GL 


a 


-0.26 


-0.25 


(-0.36,-0.16) 


(-0.43,-0.06) 


Q + /3 


0.55 


0.54 


(0.47,0.63) 


(0.40,0.71) 


se 


-2.17 


-2.17 


(-2.24, -2.10) 


(-2.31,-2.03) 




13.60 


12.96 


(9.52, 17.81) 


(7.60, 23.29) 



Table 2. Results for the fit to the full lens sample. 



Markov Chain process, we thin the chain by taking 1 point 
each 25. We find that the convergence test is well passed 
for a chain containing 10^ points which reduces to 2801 af- 
ter the burn in cut and the thinning. Such a large sample 
is more than sufficient to sample the four dimensional pa- 
rameter space. Note that we have also tried to change the 
burn in cut and thinning thus checking that the results are 
unaffected by these (somewhat arbitrary) choices. 

6.1. Constraints on the model parameters 

The results obtained fitting our model to the full SLAGS 
sample are summarized in Table 2 which reports, for each 
parameter, mean and median values and 68 and 95% confi- 
dence limits. First, we note that the best fit model param- 
eters turns out to be : 

(a, a + P,SE,tJ-vtr) = (-0.34, 0.52, -2.16, 12.22) 

giving x^/d.o.f. = 63.34/38 = 1.66. From a statisti- 
cal point of view, such a reduced value could seem 
too large. Actually, a large ^ jd.o.j. may be due to a 
model failure or to an intrinsic scatter which has not been 
taken into account. Such a second possibility is indeed the 
more realistic considering that, having stacked all the lens 
together notwithstanding their different photometric and 
mass properties, a certain degree of scatter is indeed ex- 
pected. Moreover, the overall quality of the fit appears to 
be quite good. In order to quantify the agreement between 
the data and the model for the best fit parameters, we first 
introduce the two following quantities : 

AMg _ M;ll^[iE) - M*tjiCE,Pbf) 

AcTap _ (Tap.obs — ^ap,th{Vbf ) 
^dyn ^dyn 

with pb/ the set of best fit parameters. It turns out that : 
{AME/eiens) = -0.12 , rms{AME/eE) - 1.09 , 

{Adap/Sdyn) = -0.04 , rms{A(Jap / E dyn) = 1-35 . 

The low values of {AM e / Si^ns) and {Auap/ Sdyn) are a clear 
evidence that the model is working well, while the larger 
rms values are expected considering that we are stacking 
galaxies with different intrinsic and environmental proper- 
ties. It is worth stressing, moreover, that both the projected 
mass within the Einstein ring and the aperture velocity dis- 
persion predicted by the best fit model is well within la 
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Table 3. Results for the fit to the binned lenses. 



from tlie observed values for most of the galaxies, while 
there is only one lentH with a theoretically predicted veloc- 
ity dispersion 3cr larger than the observed value. We there- 
fore safely conclude that the model have successfully fitted 
the data for the lens sample we are considering. 

Having established the validity of our model, let us now 
discuss the constraints on its parameters. As a first lesson, 
we note that models with a > are definitely excluded 
contrary to the discussion in Sect. 2 where we have shown 
that negative a values lead to a global M/L ratio diverg- 
ing for 77 — > 0. However, such an unphysical behaviour is 
actually never approached since, as we have seen in Sect. 5, 
our phenomenological ansatz do not fit realistic models for 
77 < 0.01 so that one has not to trust its extrapolation in 
the very inner regions. On the other hand, the analysis in 
Sect. 5 have shown that models with a < are needed to 
fit the global M/L ratios of different dark halo models. It is 
therefore reassuring that the likelihood analysis indeed pre- 
fer this kind of models thus leading further support to our 
phenomenological approach. As a further remark, we also 
note that the values of a and a + (3 are typical of cored dark 
halo density profiles (such as the NIS and Burkert models). 
However, the very small a values obtained when consid- 
ering the NFW model are also partly due to the adopted 
mass - concentration relation so that it is likely that chang- 
ing the c- Myir relation may lead to larger and still negative 
a values. Exploring this issue is outside our aims here, but 
we argue that cored profiles are preferred on the basis of 
our phenomenological approach. 

A further success of our model is represented by the 
constraints on s^, i.e. the logarithmic slope of the total 
density profile at the scaled Einstein radius, i.e. se = 
d log p/d log 77(77 = Re / Ref f)- Previous analysis in litera- 
ture have typically fitted the observed Einstein radius and 
velocity dispersion using a parameterized density profile. 
For instance, Koopmans & Treu (2003a, b) have fitted the 
data on 0047-281 and B1608-)-656 assuming p oc and 
finding 7 = -1.90lg;^^ and 7 = -2.03 ± 0.14, respectively. 
Using the same methodology, but averaging over a sample 
of 15 SLAGS lenses, Koopmans et al. (2006) have found 



again 7 



-2.01: 



0.03- 



As yet said before and also stressed 



This is the lens identified as SDSS J234111. 6-^000019. In 
all the fits we will consider, this is, indeed, the lens showing 
the larger normalized residuals for both the projected mass and 
the velocity dispersion so that it is possible that some peculiar 
feature is at work. However, in order to not reduce the number 
of constraints without a definite motivation, we prefer to retain 
this object in the sample even if we caution the reader againt 
overrating its large residuals. 



by the same authors, these may be considered as constraints 
on the logarithmic slope at the Einstein radius since this 
is the range mainly probed by the data. It is therefore re- 
assuring that, using a similar data analysis but a radically 
different parametrization, we get compatible constraints on 
Se thus strengthening our results. 

Finally, we consider the ^.mr parameter. It is easy to 
transform the constraints on this quantity into one on the 
DM mass fraction at the virial radius thus getting : 

{fvir} — ifvir^r^ied 0.92 , 

68% CL : (0.89,0.94) , 95% CL : (0.87,0.96) . 

Using f„r = l-espfb with /b ~ 0.17, we get {esf) ^ 0.47, 
i.e. only ~ 50% of the baryonic mass is converted into 
stars. This is consistent with semianalytic galaxy forma- 
tion models (jBenson et al. 200"0)) which typically use val- 
ues in the range (0.2, 1.0). From the observational point 
of view, the situation is much more complicated. On the 
one hand, studies matching the luminosity and mass func- 
tions (Marinoni fc Hudson 20021 |Guzik fc Seljak 2002P or 
M/L gradients ( Napolitano et al. 2005[ ) of local galaxies 
suggest esF ~ 0.25, while stellar masses in the SDSS survey 
(Padma nabhan et al. 2004) imply esp — 0.2 — 0.8. Given 
such weak observational constraints, we can only argue that 
our estimate of p,yir and hence esp is fully realistic. It is 
worth noting, however, that both observations and theo- 
retical modelling point towards a U - shaped trend of esp 
with the mass (Benson ct al. 2000, Napohtano et al. 2005] 
|Conroy fc Wec&sler 2008D . One could therefore argue that 
our approach based on a universal esp is unmotivated, but 
the small mass range probed by the lenses and the good 
results obtained let us conclude that this choice does not 
bias anyway the constraints on the quantities of interest. 

6.2. Binning the lenses 

A simple analysis of the best fit residual to the full sample 
shows that there is actually no correlation with cither the 
lens redshift, the total luminosity and stellar mass. Such 
a result may argue in favour of our assumption about the 
universality of the parameters {a, a + /S, se, Pvir)- However, 
this can also be a consequence of an erroneous estimation of 
the errors or of insufficient statistics. To further explore this 
issue, we therefore divide the lenses in four almost equally 
populated bins (denoted as Bl, B2, B3, B4) according to 
their absolute V magnitude and run our MGMG algorithm 
using chains with 2 x 10^ points reducing to 5601 after 
burn in cut and thinning. Note that, because we use now 
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Fig. 1. Variation of the model parameters with the median log Lv of each bin. The central dot denotes the median values, 
while the vertical bars span the 68% CL range. We overplot the constraints from the fit to the full sample as horizontal 
lines with the short dashed one referring to the median value and the long dashed ones delimiting the 68% CL range. 
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Table 4. Marginalized constraints on the model parameters for different bins. 



a smaller number of constraints (namely, 6 — 8 vs 38), we 
have to run longer chains in order to achieve the same con- 
vergence as with the full sample. We summarise in Table 3 
the best fit values and statistics of residuals for the fits to 
the binned samples, reporting in the second column the me- 
dian value and the full range of My in each bin, while the 
two numbers in columns 4 and 5 refer to the mean and rms 
values of the lensing and dynamics residuals, respectively. 
Table 4 reports the constraints on the model parameters 
after marginalization using the notation x^y\ to mean 
that X is the median value, while 68% and 95% CL read 
{x — yi,x + zi) and (a; — y2jX + Z2), respectively. These 
results allow us to make some interesting considerations. 

As a first remark, we note that the fit is still success- 
ful as witnessed by both the low normalized residuals and 



reduced values^. Actually, one could wonder why the re- 
duced values for bin Bl and B4 are smaller than unity 
thus arguing for a possible overestimate of the errors. While 
this is possible since our estimate is based on a naive prop- 
agation not taking into account correlations among the un- 
certainties on photometric quantities, we nevertheless note 
that the /d.o.f. values for bins B2 and B3 do not present 
such a problem. Notwithstanding the solution of this ambi- 
guity, we however stress that the constraints on the model 
parameters are not biased. 

It is interesting to compare the results on the marginal- 
ized model parameters with the corresponding ones for the 
full sample. To this aim, we look at Fig. [1] where we have 

^ The high ^ jd.o.j. for B2 is due to to the peculiar lens SDSS 
J234111. 6-1-000019 so remember the caveat yet quoted before. 
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plotted the median and 68% CL estimate of each param- 
eter as a function of the median \ogLv of the bin and 
overplotted the results from the fit to the full sample. As 
it is clear, it is not possible to infer any statistically reli- 
able trend for a which is indeed constant within the errors 
and compatible with the result from the fit to all lenses. 
A similar discussion may also be done for a -\- (3, and 
se, which due to large uncertainties cannot allow to deter- 
mine strong trends with luminosity. However, these figures 
seem to suggest interesting trends as a function of logLy, 
which we will better explore in a future paper enlarging 
the sample. Our results give issues about a dependence of 
the parameters above on luminosity and stellar mass, with 
brighter galaxies having on average higher values ior a + (3 
and fjLmr and a profile more similar to the isothermal sphere. 
Finally, an increase of with luminosity and stellar mass 
implies that the star formation efficiency esp decreases 
with these quantities. Indeed, this is consistent with the 
observational and theoretical studies ([Benson et al. 20001 
[Napolitano et al. 2005[ [Conroy fc Wechsler 20081 Tortora 
etal. 2009a, Napolitano et al. 2009) showing that more mas- 
sive galaxies are less able to convert gas into stars. 

As a general remark, it is worth noting that all the 
eventual trends seen in Fig.[l] become statistically mean- 
ingless should we have plotted the 95 instead of the 68% 
CL range. In contrast, one could equally well conclude that, 
to a first reasonable approximation, the model parameters 
(a, a + /3, Sb, Muir) do not depend on the luminosity. As 
we will show later, the derived quantities (such as the DM 
mass fractions and the scalelength parameter log 770) are 
consistent within the bins and the full sample. Motivated 
by these considerations, we therefore hereafter refer to the 
constraints from the fit to the full lens sample. We, however, 
plan to further explore this issue with a larger lens sample 
covering a wider My range in order to both improve the 
statistics and increase the number of points to infer any 
possible trend with the luminosity. 

6.3. The DM mass fraction 

Although the model we are testing has been introduced as 
an effective one to describe the galaxy, it is nevertheless an 
ideal tool to derive model independent constraints on the 
DM content of ETGs. Indeed, Eq. (|48| shows that constrain- 
ing T(r) is the same as constraining the spherical DM mass 
fraction foMir). In order to quantify these results, we sum- 
marize in Table 5, the median and 68 and 95% CL ranges 
for the projected mass fraction fprojiRs) at the Einstein 
radius and for the spherial DM mass fraction at both the 
Einstein and effective radii. To this end, we compute these 
quantities for each point of the chain obtained from the fit 
to the full lens sample after burn in cut and thinning. 

Before discussing the constraints we obtain, it is worth 
wondering whether they depend on the choice to fit the full 
lens sample or the binned subsamples. To this end, we plot 
in Fig. [2] the median values (with their 68% CL ranges) for 
the quantities of interest. It is worth noting that the errors 
bars are actually underestimated since they do not take 
into account the uncertainties on the stellar M/L ratio and 
the photometric parameters {Ie,Reff)- Indeed, the typical 
uncertainty increases in this case up to ~ 0.1 so that fully 
dominates the error induced by the model parameters. We 
have, however, adopted this choice in order to improve the 
reliability of the figures. With this caveat in mind. Fig. [2] 



shows that there is essentially a quite good agreement be- 
tween the two estimates. A direct filP^ (without taking er- 
rors into account) indeed gives : 



J proj \-^E } 



0.04-H0.93/i," (i?B) {a. 



fEMiRE) = 0.08 



: 0.02) 
0.04) 



fElliReff) - -0.003 + 0.99 f^'l^iReff) [ar 



0.03) 



Considering that the typical uncertainties are much larger 
than the offsets and that the slopes of these relations are 
close to 1, we hereafter only considers the results from the 
fit to the full sample as our estimates of the DM mass frac- 
tion^3 summarized in Table 5. 

We can further explore the viability of our results com- 
paring our estimated /z?A/(-Re//) with previous estimates 
in literature. Actually, a detailed comparison is not pos- 
sible because our systems are at an intermediate redshift 
(z ~ 0.1 — 0.5), while dynamical and X-ray studies typi- 
cally refer to nearby galaxies. Moreover, such estimates are 
model dependent relying on a given parametrization of the 
dark halo, while our constraints come from a phenomeno- 
logical and hence model independent analysis. Since it is 
not possible to take quantitatively into account these differ- 
ences, we caution the reader to not overrate any agreement 
or disagreement among our own and literature results. 

As a first preliminary step, we briefiy review some of 
the previous determinations. An accurate quantitative es- 
timate of foAiiReff) has been obtained by Gerhard et al. 
(2001) from a dynamical analysis of the line profiles of 21 
mostly luminous, slowly rotating and nearly round ETGs 
with well measured velocity dispersion and circular velocity 
profiles. According to their preferred minimum halo mod- 
els, the DM mass fraction is foAiiRef f) — 10 — 40% at 
the effective radius, finding that their central dynamical 
M/L well agree with stellar M/L estimated from synthetic 
models assuming a Salpeter IMF. Using orbit superposi- 
tion method, Thomas et al. (2005) managed to reproduce 
in detail the observed kinematics of NGC 4807 finding 
fouiReff) — 0.21±0.14 taking for a guess from popu- 
lations synthesis models with a Kroupa IMF. A larger sam- 
ple of 24 ETGs with dynamics within ^ 1 — 2Reff recon- 
structed from integral field spectroscopy was considered by 
Cappellari et al. (2006). Both two integrals analytic mod- 
els and three integral orbit superposition agree each other 
giving for /e a median value ~ 0.3 using a Kroupa IMF as 
input to the estimate of T*. Analyzing a sample of local 
galaxies, Tortora et al. (2009a) found that luminous galax- 
ies have typical values of JoMiRef}) ^ 0.3 — 0.4, using a 
Salpeter IMF, corresponding to ^ 0.6 — 0.7 for a Chabrier 
IMF. Finally, stacking the same 22 lenses we are using, 
Gavazzi et al. (2007) used a combined weak and strong 



Both in Fig. [5] and in the fits presented in this section, we 
first correct for the small asymmetry in the errors following the 
prescriptions given in D' Agostini (2004). Denoting with x the 
central value and with A± the negative and positive error so 
that the 68% GL reads (x — A_ , x + A+) , the corrected value is 
Xcorr = X ~>r (A+ — A_) with a symmetric uncertainty Acorr = 
(A+-f A_)/2. Note that, since the 68% GL ranges in Table 5 are 
indeed quite symmetric, these corrections are actually negligible 
for most of the lenses in our sample. 

The complete set of constraints, including the ones from the 
fit to the binned samples, are available on request to the authors. 
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Name fprojiRs) /DM(-Rg) fpMjRcff) VtT 

J002907.8-005550 mi^tlT, t','^, 0.564lO:gl^ 0.6021^1^9 tg.Si? 0.47t«;f, Ig;^,^ 

J015758.9-005626 0.660«:°1« 0.601«:«1^ irosl 0.6231^1? t^^, O-^&tVA life 

J021652.5-081345 0.715l»:«- t^olt 0.687«:«}? 0.7141^^6 ^ro'a^ 0.62t«:?° t";- 

J025245.2+00358 0.666«:«1^ l»:°f, 0.621«:«1J irofe 0.6511^1^ ^rofs 0.55l°:l« t'ofe 

Tnqqm9 1 nn9n'=;9 n r'^i +o"2o +0 037 rro+0.017 +0.032 „ ^„„+o.oi7 +0.032 rn+0.15 +0.29 

JUOOUlZ.l-UUZUOZ U.D01_Q (J21 -0.044 "-^"^"-O.OIS -0.036 ' "-0.018 -0.037 '-'•^'^-0.13 -0.23 

J072805.0+383526 0.563^1^1 t^!^ 0.474« °?| 0.4821^1? ti°oll 0.46t«- tg;^,? 

J080858.8+470639 0.609t°:°}^ 1°:°«3^ 0.543t° «» IS'.gal 0.5701°:°!? ti°ofr OASt^H 

J090315.2+411609 0.6781°;°- l°o:"of7 0.615l°:°l^ l°o-°f, 0.627l°:°f9 irol 0.6ll°:i« l°:rs 

J091205.3+002901 0.730t°:°}l 1° °^ 0.714«-°1« 1°;°|I 0.7321°;°;? iro'3! 0.611°;^^ 1°:« 

J095944.1+041017 0.5751°;°^? l°:°f, 0.5231°;°^ t^ot 0.5661°;°^ iro55 0.42lg:i? 1°;^,^ 

J102332.3+42300 0.5461°;°^° 1°:°1I 0.454l°:°;« 0.4651°;°;^ IS'.^ofs 0.441°;;; l°:f, 

J110308.2+532228 0.7051°°?? l"o°o|^ 0.692l°°^« l^.gf^ 0.707lg:°f9 IS.SS 0.54l°:;«3 t'.f, 

J120540.4+491029 0.622l°°;« l°:°f3 0.5731°°;^ l°°of, 0.610lg:°;«B tS.Sfo 0.48l°;;^3 

J125028.3+052349 0.5921°;°;^ l^o.Sfo 0.527l°;g;? tSfai 0.5591° °;«b irofo 0.46l°;;? 1°;^,? 

J140228.1+632133 0.6541°;°;^ l^ors 0.6051°;°;^ t'.fi, 0.6361°;°;? irofo 0.53l°;;i l":?^ 

J142015.9+601915 0.4891°;°^ IgiSfs 0.437i°;«^ ISlfo?* 0.4881°;°^^ iro's^ 0.34l°;°«9 Ig;!! 

J162746.5-00535 0.6661°;°;^ l° °i 0.6261°;°;^ t'df. 0.6591°;°;? iro!^ 0.54l°:;« 1°;^,^ 

J163028.2+452036 0.6391°;°^^ tS.Sfs 0.568l°;°;i IS.Sfe 0.5791°;°;^ l^o.S'a? 0.55l°;;i 1°;^,^ 

J230053.2+002238 0.6601°;°;^ l^o'.Sfa 0.615l°;g;i IS.gfs 0.6471°;°;? l^o.Sfo 0.54l°;;^ l°;fe 

J230321.7+14221 0.6461°;°;^ t',^^ 0.6081°;°!^ irof2 0.643l°;°o;«9 0.50l°;;^ 1°;^^ 

J234111.6+000019 0.659l°;°;i t','^, 0.6271°;°!^ IS^ofs 0.6591°;°;^ irofs 0.5ll°;;« l°;fe 



Table 5. Median values and 68 and 95% CL ranges for the DM mass fractions using the chain for the fit to the full lens 
sample. 



lensing analysis leaving the stellar M/L as a free parame- 
ter of the fit, finding out foMiReff) = 0.27±0.04. These 
authors estimate an average value for stellar M/L of 3.1 
(using h — 0.7), which is 1.5 times the median value we 
obtain using a Chabrier IMF. Renormalizing all the esti- 
mates to the Chabrier IMF, the estimates for /da/ (.Re//) 
then span the range (0.30, 0.70). Looking at Table 5, one 
can see that the values for foMiReff) span a narrow range 
with foMiReff) — 0.61 ± 0.07 averaging over the 21 lenses 
in the sample. Such a value is in agreement with local esti- 
mates even if it is worth noting that most of the lenses have 
values of foM^Ref f) at the upper end of the local range. 
Indeed, our estimates are on average consistent with the 
local ones obtaibed in Tortora et al. (2009a) when using 
the same Chabrier IMF adopted here. A more careful and 
homogeneous (with respect to the local galaxies) sample 



is needed to analyze the possible evolution of DM fraction 
with redshift which is not the case here. 

It has been recently reported in literature (Benson et 
al. 2000, Marinoni & Hudson 2002, Graham & Guzman 
2003, van den Bosch et al. 2003a, b, Napolitano et al. 2005, 
Tortora et al. 2009a, b) that the DM content correlates 
with the total luminosity, with more luminous (and mas- 
sive galaxies) having higher values of fomiRef f)- This is 
what we indeed observe in Fig. [3] where the projected and 
spherical DM mass fractions (at both the Einstein and ef- 
fective radii) are found to correlate with logLy, with this 
latter quantity estimated at the lens redshift. A direct fit, 
not taking into account the errors, gives : 

log fproj{RE) = -2.04+ 0.17 log Lv. {(J„ns = 0.03) , 

log/DA/(i?£) = -2.66 + 0.221ogLy (a™, = 0.04) , 
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Fig. 2. Comparison between different DM content estimators as constrained from the fit to the fuU or binned samples. 
Error bars refer to the 68% CL ranges. The sohd line superimposed on the data is the diagonal f{x) = x. 



log foMiReff ) = -2.41 + 0.201ogiy (ct™, = 0.04) 



which, for our model, reduces to : 



Note that the low scatter is actually far smaller than the 
typical uncertainty due to the errors on logLy and T*, be- 
ing of order 0.1. The observed correlation with the luminos- 
ity suggests that a similar correlations should be observed 
with the total stellar mass. We indeed find : 

\ogfproj{RE) = -1.64 + O.lSlogM, {ar,ns = 0.03) , 
log f DM (Re) = -1.98 + 0.15 log Af, (a,™, = 0.05) , 



logfoAiiReff) - -1.74 + 0.15 log Af, {ar 



0.05) 



where M^, denotes here the total stellar mass. Although 
encouraging, such correlations should be taken with some 
caution because of the large measurement errors on the 
mass fractions and the relatively narrow luminosity and 
mass range probed (spanning roughly only one order of 
magnitude) . It is therefore mandatory to both increase the 
sample and reduce the errors on the estimated stellar M/L 
ratio before trying any detailed comparisons with available 
theoretical scenarios. 

A different way of characterizing the galaxy DM content 
has been proposed by Napolitano et al. (2005) introducing 
the empirical DM gradient : 
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We estimate this quantity for the lenses in our sample us- 
ing the constraints from both the fit to the full lens sample 
and to the binned samples setting (rjimriout) = (0.5,4.0) 
as in Napolitano et al. (2005). A comparison between the 
two is shown in the lower right panel in Fig.O Although 
the two set of constraints are in clear disagreement if we 
consider only the median values (with a direct fit giving 
V;T(_Bzn) cx 0.83V;T(An) with a significant scatter), the 
error bars are so large that we can neglect this discrep- 
ancy and rely on the results of the fit to the full sample. 
The lower right panel in Fig. [3] then shows that this quan- 
tity is found to correlate with the luminosity. It is worth 
mentioning this trend with luminosity will be steeper if 
considering the fit binning the lenses, in agreement with 
a star formation efficiency decreasing with luminosity and 
mass ( [Napolitano et al. 2005^ . Using, however, the esti- 
mates from the fit to the full lens sample, we find : 

log V;T = -2.86 + 0.23 logLy (o-rms = 0.04) , 

logViT = -2.53 + 0.20 log Af^ (cr^™, = 0.04) . 

Fig. 6 of Napolitano et al. (2005) plots the estimated V;T 
values of a sample of local ETGs vs the total stellar mass, 
with 9.5 < log A/^ < 12. The empirical M/L gradient seems 
to be a gently falling function of log Af* over the range 9.5 < 
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log Mi, < 10.5 to become then an increasing function with a 
slope depending on the mass range considered. Indeed, over 
the mass range probed by our sample, the increase is quite 
gently (if not constant) so that our result V/T oc Af^-^° is 
in qualitative good agreement with what may be inferred 
by the local sample. It is nevertheless worth stressing that 
our V/T values are consistent with those in Napolitano et 
al. (2005). Needless to say, a larger statistics is mandatory 
in oder to confirm these results and compare them with the 
trend predicted on the basis of dark halo models. 

6.4. The impact of the IMF choice 

All the results discussed insofar rely on the initial choice of 
the stellar IMF. Adopting a different IMF, indeed, changes 
the estimated stellar M/L ratios and hence the values 
of Myir and then of Te// through Eq. (l49)) . It is therefore 
worth wondering how the results change with the IMF. To 
this aim, we rescale all the stellar M/L ratios by a factor 
1.8, that is the same as assuming a Salpeter rather than 
a Chabrier IMF. Then we repeat the fit to the full sample 
with these rescaled T* and compare the best fit parameters 
and the marginalized constraints with those in Table 2. 

Not surprisingly, we find a very good agreement for 
{a, a + I3,se), while the only parameter changed is fivir- 
This can be easily understood considering that, for a given 
lens, rescaling T* is the same as rescaling Tg// leaving un- 
altered the shape of the global M/L ratio T(r). In order to 
still fit the data, we have therefore to retain the same func- 
tional expression and hence the same values of the slope pa- 



rameters {a,a + /3, se), but scale ^eff of the same amount 
as the stellar M/L ratios. Since the luminosity is the same 
and the virial mass has to remain unchanged in order to 
give the same values of the projected mass and velocity 
dispersion, the virial DM mass fraction has to be rescaled 
according to ,: 

= 1 - T/Lt/M„. = 1 - (1 - /,^,,)(T//T,c') 

and for the mass scaling parameter we get : 

tJ.gr VT/Xt/ VT/Lt/ T/ 

in agreement with what we find within the errors. 

As a consequence, all the results obtained on the cor- 
relations between the DM mass fractions and the stellar 
luminosity and mass are still valid. Indeed, to take account 
of the change in the IMF, we have simply to rescale the /da/ 
values using /^^j — 1 — 1.8(1 — f£j^,f), while log Af* has to be 
increased by a constant factor ~ 0.25. Both these changes 
obviously rescales the corresponding zeropoints and make 
the slope analyzed steeper than those obtained using a 
Chabrier IMF strongly pointing to an increasing DM con- 
tent as a function of luminosity and stellar mass. 

7. Conclusions 

According to the estimated baryon and total matter den- 
sity parameter predicted by the cosmological concordance 
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model, DM should represent more than 80% of the total 
matter budget. As such, it is supposed to be ubiquitous 
representing therefore a significant component of any galac- 
tic structure. Notwithstanding this reasonable hypothesis, 
there is still a strong debate regarding the DM content of 
ETGs with a variety of dynamical, lensing and X - ray stud- 
ies leading to contrasting results. The lack of a reliable mass 
tracer in the ETGs outer regions and modeling uncertain- 
ties make the problem of DM in ETGs a difficult one to 
conclusively address. In order to reduce any possible sys- 
tematic bias due to the choice of the dark halo profile, a 
phenomenological approach has been proposed in Tortora 
et al. (2007) relying on the use of a versatile analytical 
expression for the global M/L ratio T(r). Such a method 
allows to smoothly interpolate between the two extreme 
cases of light traces mass and dark halo dominance mim- 
icking well the main dynamical properties of a large class 
of intermediate models. As a next mandatory step, we have 
here investigated the viability of the model by contrasting it 
against the measured Einstein angles and aperture velocity 
dispersions of a sample of 21 lenses observed by the SLAGS 
survey. A Bayesian statistical analysis makes it possible to 
obtain the main results sketched below. 

1. The effective galaxy model obtained by combining 
our ansatz for T(r) with the deprojection of the 
Sersic profile provided by the Prugniel & Simien 
(1997) approximate expression is able to reproduce 
the observations under the assumptions that the slope 
parameters {a, a -\- /3), the logarithmic slope of the 
total density profile at the scaled Einstein radius, 
se = d\np/d\nr]{r] — Re / Ref f): and the virial mass 
ratio, = Myir/T^ Lt, are universal quantities, 
while the model scalelelength 770 = ro/i?e// and global 
M/L ratio at the effective radius Tg// may change on 
a case -by -case basis. 

2. According to the marginalized constraints from the 
fit to the full lens sample, the slope parameters a 
and a -\- (3 are in agreement with what is expected in 
order the model mimics well cored halo profiles, while 
the total density profile is locally approximated by an 
isothermal model as yet found by previous studies in 
literature. We therefore argue that previous finding 
that an isothermal model fits the lensing and dynamics 
data is likely a consequence of having forced the total 
density profile to have a constant slope or of the data 
covering a too limited radial range to detect the actual 
change in the logarithimic slope. 

3. Binning the lenses according to their absolute V mag- 
nitude does not improve the fit and, on the contrary, 
gives constraints on both the model parameters and 
DM mass fractions that are in agreement (within 68 
or 95% CL ranges) with those coming out from the fit 
to the full lens sample. However, the small statistics 
and the limited luminosity range probed does not 
allow us to draw any definitive conclusion on the uni- 
versality of the proposed phenomenological T(r) model. 

4. Using stellar population synthesis models and a 
Chabrier IMF to estimate the stellar M/L ratios, 
we can constrain both the projected and spherical 
DM mass fractions at the Einstein and effective radii. 



Moreover, we have also estimated the empirical M/L 
gradients proposed in Napolitano et al. (2005) as a 
different way to quantify the DM content in ETGs. 

5. There is a clear trend of /dm (-Re//) with both logLy 
and log Mi, showing that the more luminous and/or 
massive a galaxy is, the larger is its DM content 
within the effective radius. Unfortunately, the ranges 
in log Lv and log probed are too small to make any 
quantitative comparison with theoretical models. 



Although quite interesting and in good agreement with 
previous findings, the above results are still affected by 
too large uncertainties so that it is worth wondering how 
they can be strengthened. On one hand, one can simply in- 
crease the lens sample. Indeed, while this work was nearly 
at its end, the SLAGS collaboration released a final cata- 
log (jBolton et al. 2008|) containing 170 lens candidates from 
which a subsample of 53 lenses with reliable determinations 
of the photometry, Einstein angle and aperture velocity 
dispersion. One can naively expect that more than dou- 
bling the number of objects will make our analysis more 
stringent making it possible to increase both the number 
of bins and the objects within a given bin. Unfortunately, 
more statistics does not automatically imply higher preci- 
sion. This is likely the case for our model since the data are 
still probing only a limited range centred on 0.5 — 0.7 Reff 
so that stacking more lenses does not increase the radial 
range probed. As a consequence, therefore, we do expect 
to narrow the constraints on the parameters mainly de- 
terming the shape of T(r) over this range (i.e., a and se), 
while we cannot anticipate whether a similar improvement 
may be achieved for a + P and and hence foAiiReff)- 
Therefore, rather than reducing the error on this quan- 
tity (that also depends on that on and the photomet- 
ric parameters le and Reff), we argue that such a worth 
to do analysis will improve the study of the foAiiReff)- 
Lv and fDAiiReff)-Mi, relations thus improving the com- 
parison with previous results and theoretical expecta- 
tions. It is worth noting, however, that the SLAGS sur- 
vey selection function is biased against low mass systems 
(|Bolton et al. 2006)) so that such an analysis will likely be 
unable to explore the full ETGs mass range which could 
lessen the comparison with galaxy formation scenarios. 
Enlarging the sample will allow the possible physical inter- 
pretations of the trend we obtain to be better investigated. 
The change of DM within effective radius as a function of 
mass (luminosity) scale is the main driver of Fundamental 
Plane (Tort ora et al. 2009bp . due to the change of effective 
radius (Shankar & Bernardi 2009) and/or star formation 
efficiency (Napolitano ct al. 2005,2009). Galaxy merging, 
close encounter, supernovae and AGN feedback could be 
the physical phenomena which give reason of these trends 
(Benson et al. 2000, Kaviraj et al. 2007, Khalatyan et al. 
2008, Schawinski ct al. 2008, Rusz kowski k Springel 2009, 
Shankar & Bernardi 2009, Tortora et al. 2009 a}. 

As a complementary approach, one can resort to a dif- 
ferent class of objects. Rather than concentrating on in- 
termediate redshift galaxies, nearby ETGs may be investi- 
gated relying on the velocity dispersion profile to probe 
the model in the inner regions, while planetary nebulae 
may extend deep in the outer regions thus constraining the 
outer asymptotic slope. The analysis of the best fit residuals 
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here obtained has shown that no correlation with the red- 
shift is present thus suggesting that the model parameters 
(a, a + (3 , s E , fJ-vir) do not evolve. Such a conclusion must 
be, however, checked by both dividing lenses in rcdshift bins 
and looking at a nearby sample. Should this result be con- 
firmed, one could then investigate what physical interplay 
between stellar and dark components is able to give rise to 
a global M/L ratio whose shape does not evolve with z. In 
future analysis we will analyze this problem. 

As a final comment, we would like to stress the virtue 
of a phenomenological approach (not necessarily our own) 
to the problem of DM in ETGs. A comparison with the 
determination of the dark energy equation of state w{z) is 
worth doing here. Since a plethora of theoretically moti- 
vated models provide different expressions for w{z), it has 
soon become clear that fitting all of them to the different 
data available would be a prohibitive task. This lead to the 
emergence of phenomenological proposals for w{z) able to 
explore a wide range of models in a single step also provid- 
ing model independent constraints to scrutinize different 
theories. Although the data on ETGs lensing and dynam- 
ics are still limited, it is nevertheless yet possible to rely on 
that positive experience to efficiently tackling the problem 
of DM in ETGs in a model independent way. 
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